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this end, conservative schemes which are also entropy conservative are 
constructed. These entropy conservative schemes enjoy second-order accuracy; 
moreover, they admit a particular interpretation within the finite-element 
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then show that conservative schemes are entropy stable if and only if they 
contain more viscosity than the mentioned above entropy conservative ones. 
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1. INTRODUCTION 


The systems of conservation laws referred to in the title are of the form 

(1.1) It u + £ [ f ^] = °> (x,t)C R d x[0,«); 

k=l k 

here f(^) = f(^)(u) = (f^) , • • • ,f^) ) T are smooth nonlinear flux mappings of 

r P 

the conservative variables u = u(x,t) = (up»»*,u N ) • Owing to the 
nonlinearity of the fluxes f^k) # solutions of (1.1) may develop singularities 
at a finite time after which one must admit weak solutions, i.e., those 
derived directly from the underlying integral conservative relations. 

Yet, such weak solutions of the conservative equations are not unique. 
Additional criteria are required in order to single out a unique physically 
relevant weak solution, the latter being identified as, roughly speaking, a 
stable limit of a vanishing viscosity mechansim. Entropy stability is then 
sought as the usual criterion to identify such vanishing viscosity 
solutions. Lax [10] has shown that entropy stability is in fact equivalent to 
a vanishing viscosity mechanism, at least in the small — in the large for 
scalar problems, e.g., [12], [8]. 

We study entropy stable approximations to such systems of conservation 
laws. Entropy stability manifests itself here in terms of a conservative cell 
entropy inequality. We note in passing that, if holding for large enough 
class of entropy functions, such cell entropy inequality is intimately related 
to both, the question of convergence toward a limit solution as well as the 
question of this limit solution being the unique physically relevant one, 
e.g., [1], [4], [16] and the references therein. 
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Starting with Von Neumann and Richtmyer [13] it has long been a common 
practice to ensure the entropy stability of conservative schemes by tuning 
their numerical viscosity. In this paper we provide a framework for designing 
entropy stable schemes by quantifying precisely how much numerical viscosity 
is to be added. As in [18] this is accomplished by means of comparison: we 
show that entropy stability is achieved if and only if there is more numerical 
viscosity than that present in certain entropy conservative schemes. To this 
end we proceed as follows. 

In Section 2 we begin by discussing the entropy variables associated with 
the conservative systems in question. As observed by Mock [12], [5], such 
systems are symmetrized w.r.t. these variables. Such symmetrization provides 
us with a natural order which then fits our goal to compare between the 
numerical viscosities of different schemes. Expressed in terms of these 
entropy variables, we then turn to construct the mentioned above entropy 
conservative schemes. The entropy conservative schemes — discussed in 
Section 3 and 4, are second-order accurate. Moreover, they are particularity 
attractive in view of their interpretation within the finite element framework 
which enables possible generalizations to various mesh configurations. In 
Section 5 — the main one in this paper — we compare between the viscosity 
coefficients of different conservative schemes given in their appropriate 
viscosity form . It is shown that conservative schemes containing more 
viscosity than that of an entropy stable scheme are also entropy stable. In 
particular, when compared with the previously discussed entropy conservative 
schemes, we arrive at a necessary and sufficient criterion for entropy 
stability. Since our entropy conservative schemes are second-order accurate, 
such a comparison can still entertain the construction of entropy stable 
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schemes which are second-order accurate ones* Finally, in order to simplify 
the notations, we keep our presentation within the one-dimensional case; the 
multidimensional extension can be worked out dimension by dimension* 


2m THE ENTROPY VARIABLES 

We begin our discussion with the one-dimensional model 

(2.1) yjr “ + [f(u)] = 0. 

We assume that the system (2.1) is equipped with a generalized 

Entropy Function ; a convex mapping U = U(u) augmented with entropy flux 
mapping F = F(u) such that the following compatibility relation holds 

(2.2a) U T A = F T . 

u u 

Here A = A(u) is the Jacobian matrix 


(2.2b) 


A(u) = f^. 


We note that the entropy functions, U(u), are exactly those whose positive 
Hessians U > 0 symmetrize the system (2.1) upon multiplication on the left 
[ 2 ] 

(2.3) U A = [U Al T . 

uu L uu J 



Mock [12], [5] has pointed out a more fundamental symmetrization of system 

(2.1), preserving both the strong as well as the weak solutions of the 
system. To this end one makes use of the entropy variables 


(2.4) 


/ k 3U , . 

v E v(u) (u) ' 


Thanks to the convexity of U(u) the mapping u + v is one-to-one. Hence we 
can make the change of variables u = u(v) which puts the system (2.1) in its 
equivalent symmetric form 


(2.5a) 


[ u( v) ] + [ g( v) ] = 0, g(v) = f (u(v) ) . 


3x 


The system (2,5a) is symmetric in the sense that the Jacobians of its temporal 
and spatial fluxes are 


(2.5b) H = H(v) = [u(v)j >0, B = B(v) = [g(v)]. 


This follows from the compatibility relation (2.2) equivalently expressed as 
(2-6) v T B(v) ® G(v) = F(u(v) ) , 


which in turn implies 

(2.7) u(v) =|^ [v T u(v) - U(u( v) ) ] , g(v) =— [v T g(v) - G(v)]. 
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Indeed, the Jacobians H(v) and B(v) in (2.5b) are the symmetric Hessians 
of the corresponding expressions inside the brackets in the right of (2.7). 

Finally, we note that in contrast to the symmetrization on the left quoted 
in (2.3), the use of the entropy variables symmetrize the system (2.1) on the 
right [17], i.e., the original Jacobian A = f is replaced here by the 

symmetric one B = 

(2.8) B = AU -1 = [AU _1 ] T . 

uu L uu J 


3. ENTROPY STABLE SCHEMES 

We consider conservative discretizations of the f orm^ ^ 


(3.1) 


dt \ (t) Ax ^ v+ V 2 V 2 ^ 


serving as consistent approximations to the system of conservation laws 


(3.2) 


ft u + -k t f(u) J = °- 


Here « v (t) denotes the approximation value along the gridline (x = vAx,t), 
Ax being the spatial mesh size, and 


(3.3a) 


f v+v 2 - ♦f < Vp+i'"’-'Vp ) 


' Both the differential and the discrete formulations will employ the same 
notations. The distinction between the two is made by the use of Greek 
indices in the discrete formulation. 
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is the Lipschitz continuous numerical flux consistent with the differential 


(3.3b) 


<|> f (w,w, ••*,») = f(w). 


To discuss entropy stability, we let (U,F) be any entropy pair 

T 

associated with the system (3,2). Multiplying by and employing (2.2) we 

conclude that under the smooth regime we have the additional conservation of 


entropy 


(3.4) 


+ 9F = o 
at + 3x Ut 


Taking into account the nonsmooth regime as well, following Lax [10] and 
Kruzkov [8] we postulate as an admissibility criterion an entropy stability 
requirement, expressed in terms of the following 

Entropy Inequality: We have, in the sense of distributions 


(3.5) 


aii . 3F . n 
7t + IE - 


Similarily, for the scheme (3.1) to be entropy stable, a discrete cell entropy 
inequality is sought [4] 


(3.6a) 


+7? [w 2 ‘ 


here F i, is the numerical entropy flux 


' v+ V 2 


(3.6b) 


F V+ l / 2 ’ VVq+l’-'-W 
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consistent with the differential one 

(3.6c) $ F (w,w,*»*,w) = F(w) . 

An approximate entropy equality of the above type was derived by Osher in 
[14, Section 3]. A slightly different version of such equality proved in the 
Appendix employs the consistent numerical flux. 

+ F( Vl>l +V 2l U u (a v ) * V u v + l ) ] Tf v + V 2 

(3.7a) 

_1/ 2 [ U u (u v )f(u v ) + 


The referred equality reads 


Lemma 3.1: [Osher] 

Let u (t) be the discrete solution of the conservative scheme (3.1) 
v — 

which is consistent with the system (3.2). Then, for any entropy pair (U,F) 
the following equality holds 


3r" ( "v (t » + stCv 2 - Ci/ 2 1 ' 


(3.7b) 


1 

2Ax 


u v „ Vi 

/ du U uut f v- 1/2 ~ f(u) ] + ^ dU U uuf f v+ l/ 2 " f(u) ^ 
u v-l u v 


Hence the scheme (3.1) is entropy stable if the integrals on the right of 
(3.7b) are shown to be nonpositive. In order to examine the entropy stability 
question, we first study the case in which these integrals vanish, i.e., we 
study entropy conservative schemes. 
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4. ENTROPY CONSERVATIVE SCHEMES 

In this section we identify particular schemes which satisfy a cell 
entropy equality. The numerical viscosity present in such entropy 
conservative schemes will then be used as the building block for the 
construction of entropy stable schemes. 

In order to carry out the program above, it will prove useful to work with 
the entropy variables rather than the usual conservative ones. Thus, 
associated with an entropy function U(u) are the entropy variables 
v = U u (u). Expressed in terms of the latter, the system of conservation laws 
considered is, see (2.5a) 

(4.1) [u(v)] + [g(v)] = 0, g(v) = f(u(v)). 

It is augmented by the corresponding entropy inequality 

(4 - 2a> ¥*^ 1 ° 

where (V,G) is the appropriate entropy pair 

(4.2b) V = V(v) = U(u(v) ) , G = G(v) = F(u(v)). 

In a similar manner, we interpret the conservative scheme (3.1) in terms of 
the appropriate entropy variables = U^u^). Thus, (4.1) is approximated 

V */ 2 1 


by 

(4.3a) 



where g , 1( is the numerical flux 

v+ V 2 


(4.3b) g . i, “ ♦ (v , , > • • * » v ), 

v+ V 2 g v-p+1 v+p 


consistent with the differential one 


(4.3c) ip (w,w, •••,w) 

O 


V*V #) E M 


• • *u( V ) • • • ) 

V 


“ g(w). 


The corresponding all entropy inequality takes the form 


(4.4a) 


3 T V <\<‘» + -S tV i / 2 - VVz'i 0 ' 


Here G^ + y is the numerical entropy flux 


(4.4b) 


W Vvi , " ,,T *oi ) 


consistent with the differential one 


(4.4c) 


^(wjW, • • • >w) = G(w) . 


Such a consistent numerical entropy flux suggested by (3* 7a) is given by 

Cv 2 ‘i [G<v v ) + G< Vi>] + T K + VilV'4 

(4.5a) 

• 7 [ v v g(v v > + VlS'Vl’l' 
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The equality (3.7) now simplifies into 

Lemma 4.1: (An approximte cell entropy equality-revisited). 

Let v v ( t ) be the discrete solution of the conservative scheme (4.3). 
then, for any entropy pair (V,G) the following equality holds 


3j v (v v ( t)) +■£ [g* +1/2 - g*_ 1/2 ] - 


(4.5b) 


r V 


1 

2Ax 


v+1 


/ dv [vv 2 “ g(v) ^ + / dv [ g v+ 1/2 ~ g(v) i 


lVi 


Our entropy conservative schemes will be determined by setting the 

’v+ V 2 


numerical flux g , to be 


(4.6a) 


*»+>&■ Cv 2 " 8 ( VV2 (5),d? 


where v v+ (5 ) denotes the segment connecting and 


(4.6b) 


v . 1/ (5) = v + £Av 1 , Av 1 , » v - v . 

v+ V 2 v v+ v 2 v+ V 2 v+1 v 


With this choice of numerical flux we have 


(4.7) 


v+1 


! V2 s v 2 v v 2 ' L 4 T i+ v 2 g< \ + 5 A v 1/2 )d « 


T * 


5=0 


It is here where we make use of the entropy variables formulation: thanks to 
the symmetry of g y , the expression on the right equals the path independent 
integral, see (2.7), 
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V V+1 T * 1 T Vv+1 T 

(4.S) / dv g v+i/ = /_ Av v+1 g(v v + 5Av v+1/ )d£ = J dv g(v) 

V s V 

and therefore each of the integrals on the RHS of (4.5b) vanishes in this 
case. Let us summarize what we have shown in 

Theorem 4.2: (Entropy conservative schemes). 

The conservative scheme 


(4.9) ^ [u(v v (t»] - -ji [g* +1/2 - g^ 1/2 ], g* +1/2 - NvijISWS 

is also entropy conservative, i.e., it satisfies the following cell entropy 
equality 


(4.10) 



= 0 . 


We close this section by noting that the scheme (4.9) besides being 
entropy conservative is also a second-order accurate one. Both the second- 
order accuracy and in particular the entropy conservation of the scheme (4.9) 
can be directly verified within the framework of a finite-element formulation 
advocated by, e.g., Mock [11] and in particular by Hughes and his 
collaborators , [7] and the references therein. To this end one considers the 
weak formulation of (4.1) 

(A. 11) / [|^- u(v) + -|^- g(v)]dxdt = 0. 

A 

If now the trial solutions v + v(x,t) and the weighting test functions 

A 

a) w(x,t) are chosen out of the typical finite-element set spanned by the 


C° "hat functions" 
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<x - x 


v-1 


v*> -< 


x - X , 
V v-1 


X . , - X 

v+1 

X , , - X 

V+1 V 


X , < X < X , 

V-1 — — V 


X < X < X . 

V — — v+1 


we are led to 


x 


1 a V+1 3H „ 

2 [\ + i - X v -J 97 t u(v v (t) ] = “ — g ( v(x ’ t} = 

Vl 


I v v ( t ) H v (x))dx, 
V 


or, after changing variables on the right 


(4.12) k [u(v v (t>] = x 


-2 


v+1 


- x 


v-1 


/ g (v v+ l/ 2 (5))d5 " ^ e(v ' 


5=0 


5=0 


V- 


l k 


(C))d5 


Observe that (4.12) also suggests the natural extension of the entropy 
conservative schemes (4.9) to non-unif orraily spaced meshes. The entropy 

A 

conservation follows, at once, by selecting the test function a> to conicide 

a 

with the trial function v, so that in view of (2.6), (4.11) yields 


AlJI AlJI 

0 = / [-- u(v) + g(v) ]dxdt = / [|y V(v) + -L- G(v) ]dxdt. 
SI SI 


5. NUMERICAL VISCOSITY AND ENTROPY STABILITY 

The essential role played by the numerical viscosity has long been 
recognized starting with Von Neumann and Richtmyer [13]. In this section we 
quantify the amount of numerical viscosity required in order to guarantee 
entropy stability. As in [18], this will be done by means of comparison . 

To start with, we consider the entropy conservative scheme (4.9). 
Integration by parts of its numerical flux formula (4.6a) yields 
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<W/ 2 - w 2 (E))d5 - <« ’ 

(5.1) 

■ { =0 (5 3? * ( VV 2 (5))d5 - 

Recalling the notation for the Jacobian g^, see (2.5b), 

(5.2a) B(v) = g^, 


the RHS of (5.1) can be rewritten as 


(5 - 2b) Cvt [«<V + * ( Vl ) 1 'i^ 0 <2s • 1)B(v v+V 2 (?))lJ?4v \H-V 2 • 

Inserting this into (4.9), our entropy conservative scheme assumes the 
viscosity form 


[u(v ..(t))] = - 


dt v+1 


1 

2 Ax 


[g(v 


v+1 


) " g<Vi ) ] + 


(5.3a) 


2 Ax ^v+ V 2 A V V 2 V 2 AV v- V 2 ^ ( 


Here 1 j is bhe numerical viscosity coefficient matrix given by 


(5.3b) 


Vv 2 ‘ {, 0 m - 


We note that the second brackets on the right of (5,3a) mimic a diffusive-like 

A 

terra Ax(Q v x ) x * Yet, though the viscosity martix 
not necessarily a positive definite one; rather, it is determined so as to 


l k 


is symmetric, it is 
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counterbalance the dispersive flux central differencing inside the first 
brackets on the right of (5.3a). 

Motivated by the discussion above, we would like to consider schemes given 
in a similar viscosity form [18] 


It t u(v v (t)) ] - “ 25 " g(v v-l } ] + 


(5.4) 


+ 2 Ax ^v+ V? A V V? ^v- V 2 A V V 2 ^ 


The matrix 0 i, on the right will be referred to as the numerical viscosity 
v+ a /2 

coefficient matrix. 

Remark : We observe that the above definition of the numerical viscosity 

coefficient depends on the specific entropy function under consideration 

2 

U(u). The special choice U(u) = u corresponds to our earlier viscosity 
definition [18] in the scalar case. 

What schemes admit a viscosity form-like (5.4)? To answer this question 
we observe that the numerical flux determined by (5.4) 


(5.5) 


vi/ 2 a Vv P+ i ),,, *Vp ) = i t g( 


V + g( Vl 


>3 



satisfies the consistency relation 


(5.6) 


V , v-p+l ,, "’Vl' w > w > , v+2-" , -Vp ) ■ g(,) - 
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The consistency relation (5.6) is slightly more stringent than the usual one 
(4.3c). It characterizes the class of essentially three-point schemes , a 
class which includes, beside the standard three-point schemes, most of the 
recently constructed second-order accurate TVD schemes, e.g., [3], [15], 
[6]. In the Appendix we show that the converse of the above implication 
holds, namely we have 

Lemma 5.1: The consistent conservative scheme (4.3) can be rewritten in 
the viscous form (5.4) if and only if it is an essentially three-point scheme. 


Granted the viscosity form (5.4), we now turn to discuss the question of 
entropy stability. We first note that there is only one degree of freedom in 
setting up the viscosity form (5.4), that is, setting up the viscosity 
coefficient Q^ + y • Motivated by Lemma 4.1, we shall examine the quantity in 
light of the integrals on the RHS of (4.5b). Thus, using (5.5) we find 


T 1 T 

J dv [%+ 1/ 2 - z^)} = - l/o + sK+i >1 


2 av v+ 1/2 


- J dv T g(v) - I Av* + 1/2 Q v+ 1/2 Av^ 1/2 . 

V 


The first two terms on the right are fixed by the differential flux; 

abbreviate their difference by K i, . Then, Lemma 4.1 tells us that 

v+ 72 


1 r * 


— V( v (t)) + — [G , i, - G i , 1 = K i,+K i, 
dt v Ax L v+ V 2 v— V2 J v~ V2 v+ ^2 


2 AV v- V2 ^v- V2 A V V2 2 A V V2 ^v+ V2 A V V2 


(5.7) 
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Hence we conclude that the cell expression on the left decreases whenever the 

T 

quadratic forms ^ V v ± I /2 V 2 V 2 * ncrease * This suggests to compare 
between different schemes in terms of their numerical viscosity coefficients 
used as our scale* 

We shall say one scheme contains more viscosity than another scheme if the 

viscosity coefficient of the first scheme, say Q^+\/ > dominates that of the 
( 2 ) 

second one, \j > i.e*, if we have 


(5.8a) 


T (2) T (1) 

V 2 V */ 2 4 V V 2 4 V l / 2 V 1/2 4 V 1/2 


An strengthened formulation of this, in terms of the natural order 


among 


symmetric matrices, reads 
(5.8b) 




In case the entropy variables were used to begin with, they would lead to 
symmetric viscosity coefficients very much the same way they led to symmetric 
Jacoians, e.g., (5.3b), and we would arrive at the natural hierarchy 


(5.8c) 


v+ V2 — V+ V2 


The last three inequalities mean the same when dealing with the scalar case. 

Equipped with the above terminology we now turn to a particular comparsion 
with the entropy conservative scheme (5.3). Thus, we consider a conservative 
scheme in its viscous form (5.4) and we let D . 1 , 


^v+ V2 ^v+ V2 ^v+ V2 


(5.9) 
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denotes its viscosity deviation from the entropy conservative one (5.3). Then 
we can write the scheme under conservation (5.4) 


S [" c V t)) ] = IS + IS tv V 2 ly v+V 2 ' Cv 2 4 vv 2 1 

+ 2 Ax t^v+ V 2 AV v+ V 2 D v- V 2 ^ w v- V 2 t 

Ax ^ V V2 ®v- V2 ^ 2 Ax ^v+ V2 A V V2 ^v— V2 A V V2 ^ * 


Multiplying the last equality by U^(u y (t)) = v^(t) on the left yields 

— V(v (t) + — v T [g* 1 - g* 1 1 

dt v' Ax v *■ ®v+ V2 ®v- V2 

(5.10) 

1 T 

2 Ax V y ^ v+ V2 A V V2 V2 AV v - V2 ^ * 


By Theorem 4.2 the second expression on the RHS of (5.10) equals the 
conservative difference 


1 T r * * n 

Ax V v^v+V 2 ®v- Vz 

(5.11) 


■ v v (t) 3t [“ < V t,) l 

‘it v(T v (t>) - S tvv 2 ‘ Cv 2 l 


Regarding the RHS of (5.10), the following identity puts it as the sum of 
familiar quadratic terms plus a conservative difference 
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1 T 

-r v [ D i / Av i , - D i , Av i , I 

2Ax v L v+ V 2 v+ V2 v ” V2 v- /2 J 

(5.12) = - - 57 -- [~ Av T i, D i 7 Av l, + Av T i, D 

2Ax L 2 v+ V2 v+ V 2 v+ V2 2 v- 72 v- 

+ 2^ t| K + Vl 1 T V V 2 4 V l/ 2 " £ t Vl 

Using the last equality we conclude the main result of this section, 
asserting 


V2 A v V2 ^ 

+ V V ^ ^v- V 2 AV V“ V 2 ^ * 


Theorem 5.2: (Entropy stability). 

A conservative scheme is entropy stable if and only if it contains more 
numerical viscosity than the entropy conservative scheme. Moreover, the 
entropy dissipates in this case at rate governed by the cell-estimate 

It v <V t » + t G v+ 1/2 ' G v- 1/2 1 

(5.13a) 

IT T 

= - -T— fiv . ii D . i< 4v . i, + Av 1 , D 1 / Av l, 1 < 0. 
4Ax L v+ V 2 v+ V 2 v+ V 2 v- 72 v- 72 v- 72 J — 

Here G^ + y is the consistent numerical entropy flux 

V v 2 - i + G( Vi» + 1 K + v v+ 1 ] t (w 2 ' 1 Vv/v v 2 1 

(5.13b) 

"7 C v v g(v v ) + Vl^Vl^* 

Proof: Inserting (5.11) and (5.12) into (5.10) we obtain 
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(5.14) 


It V(v v (t>) + s [°v+ 1/2 - G v- v 2 ] 

AAx ^ iV \>+ V2 D v+ V2 iV v+ V2 + iV v- V2 D ' J ~ V2 iV v- V2 ^ 


where G^ + y is the consistent numerical entropy flux given in (5.13b). Now, 
if the scheme considered contains more viscosity than the entropy conservative 
scheme, then by (5.9) the RHS of (5.14) is nonpositive and entropy stability 
follows. Conversely, assume that our scheme is entropy stable, i.e., that it 
satisfies a consistent all entropy inequality of the form 


+ i;[Vy 2 - Vvz 1 ^ 0 - 


Subtracting this from (5.14) we get after multiplication by Ax 


(5.15a) H 1 , — H 1 . > - - 7 - [Av T , 1 , D 1 . Av 1 , + Av 1 1 , D 1 , Av ]. 

v+ V2 v- 72“ 4 L v+ 72 v+ 72 v+ 72 v- 72 v- i2 v_ >2 J 


Here H 1 . = H(v v . ) stands for the entropy flux difference 

v+ 72 v-q+ 1 ’ ’ v+q 

G v+ 1 - G y+ y ; by Lemma 5.1 the entropy flux G y+ is essentially a three- 
point one and a similar form of G y+ y implies that H y+ y 
essentially three-point consistency relation, see (5.6), 


satisfies the 


(5.15b) 


H( v ,,,♦•• ,v , ,w,w, v . „ , •••,▼. ) = 0 
v-q+ 1 ’ v -1 v+ 2 ’ ’ v+q 


Choosing v v _^ = v v in (5.15a) we obtain, in view of (5.15b) 


(5.16). 


^v+ V 2 ~ 4 Av v+ V 2 °v+ V 2 A V V 2 * 
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Similarly, taking = v^ +1 in (5.15a) yields, in view of (5.15b) 

1 T 

H i , > - t 4v ], D 1 / &v i, . 
v- V 2 “ 4 v- V2 v- V2 v- V2 

Adding (5.16)^ together with (5.17)^^ implies 

® ^v+ V2 ^v+ V2 — 2 A V V2 ^v+ V2 A V V2 

or, according to (5.9), 

(5.18) 4 V l /2 V l/ 2 4 V l/ 2 < 4 ’w 1/2 V 1/2 4 V V 2 ' 

Thus the scheme considered, (5.4), contains more viscosity than the entropy 
conservative one, (5.3), as asserted. 

We recall that the entropy conservative scheme (4.9) is a second-order 
accurate one. Hence, Theorem 5.2 allows us, in particular, to tune additional 
numerical viscosity so that we retain both the entropy stability and second- 
order accuracy [15], [16]. 

Finally, using arguments similar to those employed in the proof of the 
last theorem, we conclude with the following extension of [18, Theorem 6.1], 
dealing with systems of conservation laws. 

Theorem 5.3: (Entropy stability by comparison). 

Conservative schemes containing more viscosity than an entropy stable 
scheme are also entropy stable. 


(5.17) v 
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The last theorem provides us with a framework for designing entropy stable 
schemes. One could start with any scheme which is already known to be entropy 
stable and then tune in additional numerical viscosity according to the above 
guidelines so as to obtain a better performance in terms of, e.g, simplicity, 
controlled variation, avoiding entropy "glitches" etc. Virtually all the 
first-order entropy stable schemes fall within that framework when compared to 
Godunov scheme — these are the E schemes discussed in [14], [18]. 
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APPENDIX 


A. AN APPROXIMATE CELL ENTROPY EQUALITY 

We consider the scheme 


(a.l) 


d 

dt 


U v (t) ’ ' A [**. l/ 2 - Wjl' 


Multiplication by U u (u v (t)) on the left gives us 


(a. 2) 


It u( "v (t)) ■ - A 1/2 ' U >v )£ v- V 2 )■ 


Adding the term — [ F^ + ^ - F^_ y ] to both sides, we find on account of 
(3.7a) 


(a. 3) 


+ A t F * +1/2 - Cv 2 l - m [( p <Vi> - F( Vi > ) 


2Ax 


tVvi> - V“v» Tf vH-v 2 + iw - WidVvi 


Integration by parts of the compatibility relation (2.2) implies that the 
first brackets on the right of (a. 3) equal 
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F(u) - U*(u)f(u) 


(a. 4) 


u 11 u 

v+1 v+1 I v+1 

= / F du - U (u)f(u) 

J u u 

11 , u . 'u , 

v-1 v— 1 v-1 


Vi ,Vi 

■ / Vu d " - U u«"> 

U v-I Vl 


u v+i 

-/ duu uu f(«), 

Vl 


The second brackets on the right of (a. 3) can be rewritten as 


u 


(a. 5) lT(u) 


v+1 


u 


f . 1/ - U T (u) 
v+ V2 u 



Inserting (a. 4) and (a. 5) into (a. 3) yields the desired approximate cell 
entropy equality, see Lemma 3.1, 


St u <“v (t)> 


_1 r * 

Ax l V+ 1/2 



1 

2Ax 




The above proof essentially follows that 
differs, however, in that it employs a 
centered at mesh midpoints. 


of Osher in [16, Section 3J. Il 

* 

numerical entropy flux, F 1 , , 

v+ 72 


B. THE VISCOSITY FORM OF ESSENTIALLY THREE-POINT SCHEMES 

We consider conservative schemes in the viscous form 


Sf f u(v v (t) l = aF fvv 2 " V v 2 ] = 2^ [rfW “ g( Vi } ] 

(b.l) 

+ 2 Ax ^v+ V 2 A v V 2 ^v- V 2 A V V 2 ^ 
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The numerical flux associated with scheme (b.l) 


(b.2) 


®v+ V 2 V V“P+1 * V V+p ^ 


- 7 [b(\) + ^Vi'l ■ T Vv 2 tVi - V J 


satisfies the essentially three-point consistency relation 


(b.3) 


K ( v +1 ,---,v v _ 1 ,w,w,v v+2 ,...,v ) = g(w). 


Conversely, consider a conservative scheme 


(b.4) 


[ u( v (t)l = - ~ U (v ) - Ip ( V , , )1 

L v J Ax Lr g v-p+1 v+p g v~p v+p-1 J 


Subject to the essentially three-point consistency relation 


Cb.5) 


+g ( Vp+l’ - ">Vl’"'"’V2-"-Vp ) - 8 ( »> 


the scheme (b.4) can be put into the viscosity form (b.l) provided a numerical 
viscosity coefficient, Q^ + y > can be found such that 


(b.6) Q ^1/ Av xl/ a g( v ) + g(v ) - 2$ (v ). 

v+ V 2 v+ V 2 v v+1 g v-p+1 v+p 


Using the consistency relation (b.5) we can rewrite the RHS of (b.6) as 
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+ 1* (•••.v .v . • • O — dj C • • • *v .v . • • O 1 

L V * v+r v+r V ’ v* v+r J 


and equality (b.6) is then fulfilled by setting the numerical viscosity 
coefficient to be 


(b.7) 


^>+ 1/2 


= ' u t V'”’’vW 2 (5)i, " )]dc 


+ Lo t V"' ,V v + V 2 ( 6 ),T v+l , " ) ] d «- 
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